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Abstract: 

We compute the next-to-leading order term in the long-distance expansion of the 
mutual information for free scalars in three space-time dimensions. The geometry 
considered is two disjoint disks separated by a distance r between their centers. No 
evidence for non-analyticity in the Renyi parameter n for the continuation n —> 1 in 
the next-to-leading order term is found. 
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1. Introduction 

Cardy [1] has presented a general framework for the mutual Renyi information of 
two disjoint compact spatial regions A and B for a d + 1 dimensional conformal field 
theory in the limit when the separation r between A and B is much greater than 
the sizes Ra and Rb■ The method involves the replica trick [2,3] for computing 
the Renyi entropies S^ = (1 — n) -1 logTrp^, and from these quantities the von 
Neumann entropy Sa = — Trp A log pa of the reduced density matrix p A , obtained as 

S A = lim . (1.1) 

71 —> 1 

This limit generally involves a non-trivial analytic continuation of the Renyi para¬ 
meter n. 

Similarly, from the Renyi mutual information 


/<")(As) = si" ) + 4” ) -s!.uB 


(1.2) 




















one can obtain its corresponding mutual information 


I(A, B ) = S(A) + S(B ) - S(A U B ), 


(1.3) 


taking the same limit. 

For a free scalar theory in (d + 1) dimensions, one has in general the leading 
term in the large r expansion 

1 . (1.4) 


For the simple case where A and B are spheres of radii Ra and Rb, Cardy found for 
3 + 1 dimensions 



n 4 - 1 
15n 3 (n — 1) ’ 


(1.5) 



( 1 . 6 ) 


for 2 + 1 dimensions 1 and a known definite integral for . 

A compound system with a cusp, consisting of two identical spheres in contact 
widely separated from two other identical spheres in contact also gives a mutual 
Renyi entropy (|1,4|) with [5] 
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where now in ( |1.4| ) Ra is the radius of one of the spheres in A and similarly for B. 
Thus, the leading term of the expansion is of universal form, with g^‘ } encoding the 
geometry of A and B. 

For d > 2, higher order terms in the large distance expansion of the mutual 
information I(A,B) requires the contribution of the stress tensor T^ u and current 
to evaluate the corresponding coefficients. For example, for d = 2, which is the 
subject of this work, contributions from T^ u and d^cj) only appear starting at the next- 
to-next to leading order ~ ( Ra ^ b ) and their contribution is of equal importance as 
that of (analogous to the one computed in Section 4). For d — 3, this additional 

1 Recently, the coefficient for any dimension was found to be , for 

arbitrary CFTs [4]. Here A is the lowest scaling dimension of the CFT operators and Aa its 
degeneracy. Notice that for free CFTs, the d dependence comes from the relation between A and d 
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terms contribute already in the next to leading order term and their contribution 
dominate for d > 3. 

Care must be taken in these calculations to ensure that spurious singularities 
do not appear in the many integrations to be carried out. Given the complexity 
of the calculation, one may be concerned with a possible non-analyticity in the 
Renyi parameter n. However, no evidence for such non-analyticity is found in the 
continuation n —» 1 of the next to leading order term in 1/r for free scalar fields in 
3 space-time dimensions (d = 2). 

For d — 2 free scalar fields, we have the long distance expansion 


/ RaRb\ 
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( 1 . 8 ) 


where the sub-leading term will be the subject of our analysis. It is found that 


I(A,B) 
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The appearance of (1 /vr 2 ) in (|1.9|) is unexpected. 


(1.9) 


2. Review of Cardy’s paper [1] 

Consider two disjoint compact objects A and B in a CFT, It was shown in [6] [7] [8] 
for 1+1 dimensions, that the Renyi mutual information given by 

/<”>(4B) = +” ) +4” ) -S2b, (2.1) 


has the following expansion 

I M (A,B) +C?(ft))C«({k j )) 

{+} 


(RaRbY 11 ^ h 

V r 2 ) 


( 2 . 2 ) 


where the are the scaling dimension of the operators This expansion was 


studied in detail in [8]. Cardy generalized the equation (|2.2[) for higher space-time 
dimensions, and it is the intent of this section to present a brief review of his argu¬ 
ments. 

The Renyi entropy for a given region X is given by 


o(n) _ 
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logTi 


(2.3) 


which in turn can be computed in terms of a path-integral on a conifold 
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Here Z is the partition function on the original space, and Z(C%) is the n th 
associated conifold. The Renyi entropies are then 

Z(C?3b)Z" 


J (n) (A,P) = 


1 — n 


log 


Z(Cf)Z(C' B 

The basic idea that leads to the expansion ( ]2.2| ) is the observation that the sewing 
operation on the entangling regions can be thought of by a distant observer as a semi¬ 
local operation, and so can be implemented through a weighted sum of a product of 
local operators. That is 
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The coefficients of the expansion can be obtained by 
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( 2 , 8 ) 


and the ratio of partition functions relevant for the evaluation of the mutual infor¬ 
mation is 


Z(C 
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(2.9) 


One can write the left-hand side of ( j2.9|) as 

pH 

p(n) — r A\JB 


p( n ) p( n ) 

J A r B 


( 2 . 10 ) 


where P^' 1 = E(C]”' , )/Z n , and expand ( ]2.5|) systematically in the large r expansion; 
so that 

H4 B) = - i[P<% + [P<">]( S) + ...}, (2.11) 

where [pA )corresponds to the order m term in the expansion of (|2.10| ), that is 
[P {n) U rsj (. R A R B /r 2 ) m . The [P<")] (o) term is equal to one, as it is obtained from 
the contribution of identity operator = 1, and does not contribute to Q2.11|) . 

If one is only interested in the mutual information, that is 


’Hi 


I(A,B) = lim I {n \A,B) 


n —>1 


( 2 . 12 ) 


one can ignore the second term in (|2 .1 1|) as [P (n )](i) is O^n — l) and therefore [pH]^ 
is 0(n — l) 2 . Also each term [P^]( m ) for m ^ 2 should go to zero at least as 0(n — l) 2 
when n goes to 1. 
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3. Leading contribution to mutual information 


Cardy [1] computed the leading long distance contribution to the mutual information 
of widely separated disks. In this section we will briefly review his calculation. The 
lowest order coefficients y are 

Cfji = C* 0 i i o — lim (x 1 x 2 ){(j)j(x 1 )(j) f (x2))An) (3.1) 

JJ X\,X2^00 ^A 

for j t - f and 



In the case of spheres we can use a conformal transformation that takes the spher¬ 
ical surface of interest (say one with radius R A ) in to an infinite plane R d_1 , and 
therefore the original conifold is transformed to C'[^ = {2-dimensional cone of open¬ 
ing angle 27m} x M d_1 . This transformation also takes points at infinity and maps 
them to points on a unit sphere so that the previous coefficients in terms of this new 
coordinates becomes 


Cjjt = 2R A (4>j(l)(j)j'(l)) c ,(n) , 
Cfi = 2 ~ 1/2 R a (: $(1) . 


(3.3) 


The Green functions for free scalar fields for d = 2, were constructed using the 
method of images and an analytic continuation in n. The final result is 


G (n) (l,0,O) = 


oo x ( 0 / 2 tt ) — 1 (-| _ 


(1 ~x) 


27r J o (1 + x)(l — x n ) 


dx. 


In terms of the new coordinates 


(0j( 1 )<M 1 )) c (;•) = tt; 


00 x \i-i' l-i 


(1 ~x) 


dx. 


(3.4) 


(3.5) 


c a 27r J o (1 + x)(l — x n ) 

These coefficients are what is needed to evaluate the ratio of the partition func 


tion (p.9|) , so the leading term is 


r 2 (n X C ti' C fr + X y’ 


(3,6) 


j¥j' 3 

which for the Renyi mutual information implies 2 [1] 
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(3.7) 


2 The symbol ~ accounts for the fact that the RHS of (3.7) is the leading large distance term for 
n close to 1. 








We can fix one index, and put an overall factor of n using the cyclic symmetry of 
the sum 


l( n \A,B) 


n 2 RaRb ( 


n -1 


71—1 


y ^ ' {^(l)^o(^))^M (0j(l)0o(l)) c W 


3 = 1 

jl2/i \ ,\ /. -2/ 


+2( : 0o(l) : )cw( : 0o(l) : )c^ n) ) ‘ 


(3.8) 


Since we are interested in the mutual information, we focus on the terms that survive 
the n —> 1 limit. Taking into account the fact that (: </>q(1) :) ~ 0(n — 1), it is clear 
that those terms do not contribute to (|3.8|) and can be safely dropped. The remaining 
sum is 

n— 1 

L -* '-'A 

3 = 1 
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r»QO /»00 


_ (1 - ( xy) n *)(1 -x)(l -y) 

4tt 2 J 0 J 0 (1 + x)(l + y)(l - xy)(l - x n )(l y n ) 
n — 1 f°° log (xy) 


4vr 2 J 0 Jo (l + a;)(l + y)(l-x 2 /) 


dxdy. 


dxdy 

(3.9) 


The second line in (|3.9|) is obtained from the first when the n —> 1 limit is 
taken. This integral can be done by a change of variables y' = y and x' = xy which 
transforms it to 


71 — 1 
47T 2 


'”<fa' log(x ' 
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(l-x')Jo (1 + y')(x' + y') 47r 2 

= (n - 1) 

6 

The mutual information can be evaluated by expanding the log in 
n —y 1 limit, that is 


r °° dx' {l ° g{x ' ))2 


(1-x') 2 


(3.10) 
and taking 
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4. Next to leading term 


In this section we apply the procedure of the previous section to evaluate the next 
to leading order term in the expansion (|2.11|) . 

In terms of the C-W, coefficients, the next to leading term is 
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Using the same argument as above we can do one of the sums explicitly due to the 
cyclic symmetry of the labels 
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71 — 1 


yi c tjoo c ?joo+ 


3 = 1 


n—l 


c tj'oo c fj'oo+ 


j¥=f 


n—l 
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/oA syB \ 


(4.2) 


We would like to evaluate each of these terms individually, considering only the 
0(n — 1) contributions, which will be the only ones which eventually will contribute 
to the mutual information. 

The coefficients appearing in ( |4.2| ) in terms of correlators of the fundamental 
fields and <3>j, are given by 3 

C}j oo = (2JU) 2 <Sj(1)*o(1)>, 

Cffoo = (2fl x ) 2 <4,(l)^.(l)*o(l)), 

Off ,„ = (2^) 2 (0 j (l)0 j ,(l)*(l)^o(l)) . (4.3) 

where %(1) = -jj : tf(l) :. 

The expression in ([T2]) contains terms that differ in the number of nested sums 
involved. The higher this number is the more complex the evaluation turns out to 
be. In the rest of the paper we will evaluate each of this terms 


4.1 Evaluation of £s'E'j=l c fjoo c ?joo 
Using Wick’s theorem the coefficient 

Cf j00 = 2 R\{-. $( 1 ) :: $( 1 ) :) 

= 2i£«: <j>]{l) :)(: ^g(l) :) + 2(^(l)0 o (l))(^(l)0 o (l))) • (4.4) 


In the n —> 1 limit, any one point function operator in the conifold of singularities 
goes to zero at least as 0{n — 1). As it will be evident later, any sum gives rise 
to an (n — 1) factor. Therefore, we can safely neglect any term which is already 
0(n — 1) or higher before doing a sum, in our case at hand we can ignore the term 

(: °j( 1) :)(: 4>l{ 1) :> in 0. 

That means that our term of interest is 
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2H 
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^(^(i)?i„(i)) 4 . 
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(4.5) 


This term can be evaluated using the same technique used in the evaluation of the 
leading contribution (|3.9| ) . This entails the use of the explicit integral representation 

3 For notational convenience, here after we use simple brackets () instead of the more accurate 

Oca • 
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of the two point, and carrying out the sum explicitly 


n— 1 

EW'jM'M 1 )} 4 

3 = 1 

1 roo roo (1 — (xyzw) n ~ 1 ) (1 — x) (1 — y) (1 — z) (1 — w)dxdydzdw 

167T 4 J 0 J 0 (1 + x)(l + y){ 1 + z)(l + w)(l — xyzw)( 1 — x n )(l — y n )( 1 — z n )( 1 — w n ) ’ 

(4.6) 


which in the n —> 1 limit evaluates to 4 




and so we conclude that 


r A r B - 4n ( n ~ x ) ( r a r b 
2 r 4 Z-/ °iioo°ijoo ^ 1 r 2 

j=l v 


(4.7) 


(4.8) 


4.2 Evaluation of CfjwCjj'X 

The second term in ( [f.2|) requires the evaluation of a double sum and therefore some 
extra care should be taken. First, we need Cyyoo 

Cfj 'oo = (2i?^) 2 —^=(0j(l)0j/(l) : </> 0 (l) :), (4-9) 

which can be rewritten as 

qfoo = 4J$-h«* j( l)^ (1 ))(: ^(l) :) + 2(^(1)^„(1))(^.(1)^„(1))) (4.10) 

using Wick’s theorem. 

Here again we neglect the term which contains the one point function (: </>d(l) :) 
for the same reason given above, and so we have 
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3p 
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E /^iA siB 
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oo 


16nR A R Bj2(Mi)Mi)nhmMi)Y- 


3 r 4 




32nR\R B 
3 r 4 


j-i 


5>,(l)<A„(l)) 2 £«y(l)'Ml)> 2 , (4.H) 


j=2 


i'=i 


where in the second line we have separated the sum in two, each one with a different 
order: j > f and j' > j, and used the symmetry j -H- j' to equate the two terms and 


4 See Appendix (|A]) for the details of the evaluation 










add them up. 


n— 1 


3~ 1 


^(^(i)4(i)> 2 x;<^(i)*(i)> 2 

j=2 j '=1 


167T 4 


n—1 

E 

1=2 ' 


dxdydzdw(l — (zw) j l )(xy)i X (1 — x)(l — y)(l — z)(l — w) 

(1 — zw)( 1 + x)(l + y)(l + z)(l + u>)(l — x n )(l — y n )(l — z n )(l — w n ) ’ 

(4.12) 


in the n —> 1 limit, we can evaluate exactly this expression as shown in Appendix 
where we find that 


n—i i- 1 / _i\ 


j= 2 j'=l 

Therefore, the net contribution from the double sum term is 


n—1 

n r A r<B 


8 n(n — 1) / RaR 


3!r 4 


oo 


o+i' 


45 


i B 


4.3 Evaluation of ^ Ej'yyy; C fj'io C fj'io 

The final term we want to evaluate involves the triple sum 


n 

4! r 4 


n—1 

E ^-yA ^yB 


w ■ 


(4.13) 


(4.14) 


(4.15) 


As before,we rewrite higher point function into lower ones using Wick’s theorem 

q*, i0 = 

= 4J$M^(l)^(l)X*(l)*o(l)> + 

(^(l)*(l))(^,..(1)4(1)) + <0j(l)0o(l)><0j'(l)0Kl))]. (416) 


Thus ( [4.15|) is 
16n ( RaRb 


n—1 


£ 3(0,(l)0j'(l)) 2 (0;(l)0o(l)) 2 + 


n—1 


- ^ 6(^o(l)0j(l))(^j(l)0*( 1 ))(0j'( 1 )0i( 1 ))(^j'( 1 )^o(l)) 


(4.17) 


after relabelling the indices to identify equal terms. The full expression preserves the 
symmetry j <(->■ j' so we can give an order to these two indices (j > j') and multiply 
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the final expression by two. To evaluate this we have to sum over all possible orders 
that respect the given order (j > j'). That means we will have three different terms 
from each sum, those are: 1) j > j' > l, 2) l > j > j' and 3) j > l > j'. To avoid 
confusion let use k instead of j'. The triple sum terms are 


p \ 2 n—l j— 1 fc—1 

r EEE Wj(i)^(i)) 2 (*(i)*(i)) 2 + mDummuDf 

' j =3 k =2 1=1 

+ (^(l)^(l)) 2 (^(l)0o(l)) 2 

+2(0o(l)0j(l))(0i(l)0K 1 ))(0fc( 1 )</ , i( 1 ))(0fc( 1 )0o(l)) 

+2(0o(l)0fc(l))(0fc(l)0j(l))(0i(l)0/(l))(^( 1 )</ ) o(l)) 

(4.18) 


Writing all propagators in terms of integrals, and carrying out the sums in a con¬ 
venient order, we find a similar looking expression to the ones we have evaluated in 
Sec 2, with the final answer 


71—1 

4 FT S C fj'lO C fj'lO = 4n ( n 



12 1 \ 
30 + 45 + 3^V ’ 


(4.19) 


For the specifics of this evaluation the reader is advised to look at Appendix (^). 

Adding all the evaluated terms, that is the terms with single, double and triple 
sums in (|T2|), we get 


71 — 1 

n , 1 v—-\ 


71—1 
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1 _ _ 1 ___ 

r'lB 1 \ A s~iA iB . 1 \ A /-iA r'lB \ 

r 4^2 7 -> WjOO'+j'jOO ' qi / ; ^jj'OO^jj’OO w / , ^jj'io^jj'io) 

3 = 1 
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n(n — 1) ^RaRb ^6^4 
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(4.20) 


The mutual information to second order in the long distance expansion parameter 
( RaRb ) j g then given by: 




( RaRb\ 

J 


1 / RaRb 
+ 3 \ A 


2 



(4.21) 


where the evaluation of the sub-leading term is the main result of this work. 


5. Discussion 

In this work we have presented a detailed evaluation of the next to leading order term 
in the mutual information of disks in a (2+1) space-time dimensional free field theory. 
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This coefficient has not been evaluated before, neither analytically nor numerically, 
so there are no other results which we could compare it with. However following 
the spirit of [4] one could wonder whether this coefficient is in anyway related to the 
subleading coefficient computed for CFT 2 in [9] (with a A associated to the free field 
theory case for 3 spacetime dimensions, that is A = 1/2). In fact, the numeric value 
of these numbers are surprisingly similar (for CFT 2 this is 8/15 ~ 0.533, while for 
(2+1) free field theory is 6/15 + 4/(37 t 2 ) ~ 0.5350 • • •). Although, we do not expect 
any kind of universality for the next to leading term 5 , it would be interesting to 
explore whether this approximate equality comes from some unknown bound that 
explains it. 

We can estimate a hypothetical critical radius which might be interpreted as a 
break down of the power series expansion, by estimating the distance r c at which the 
leading and next to leading order terms are of the same order of magnitude. That is 


f—) 

v RaRb J 



Assuming Ra = Rb = R, we get: 


r c ~ R 




1.267/2. 


(5.1) 


(5.2) 


This is clearly beyond the validity of the expansion, since it conflicts with the ge¬ 
ometric set up. We interpret this as a justification of the expansion (|1.8|) for large 
r » R. 

One of the motivations of this study was to explore the possibility of non- 
analyticity in the Renyi parameter n for the sub-leading term in the continuation 
n —> 1. There is no indication of non-analyticity in this continuation to the order in 
the large distance expansion we considered. 

One can also consider the holographic expansion of the mutual information of 
two regions A and B on the boundary surface. The leading term at 0(N 2 ) in the 1/A 
or Gn expansion vanishes, by the Ryu-Takayanagi argument [10]. In this context 
the leading large distance term is 0(1) rs_/ O(G 0 N ) [6] [11] [12], and a phase-transition 
between bulk surfaces of different topologies occurs [6]. These results present a chal¬ 
lenge for the holographic calculation of mutual information to understand comparable 
terms in a large r expansion 6 . 


5 As it occurs for the leading coefficient in this expansion. This is explained by the universality 
of the 2 point functions in CFT’s [4]. 

6 Notice that an exact match between the large distance holographic mutual information and its 
CFT dual was found by explicit calculation in [4], so it would be interesting to explore the next to 
leading terms in both the CFT and its bulk dual. 
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A. Calculation of single sum term 


We need to evaluate the expression (|4.6|), this is 


n— 1 

I>j( i)Mi» 4 

3 = 1 

1 f°° f°° (1 — (xyzw) 71 - 1 )^ — x)(l — y)(l — z)(l — w)dxdydzdw 

167T 4 J 0 J 0 (1 + x)(l + y)( 1 + z)(l + w)(l — xyzw)( 1 — x n )(l — y n )( 1 — z n )( 1 — w n ) 

(A.l) 


In the n —> 1 limit, this is: 

(n — 1) f°° f°° log (xyzw)dxdydzdw 

167T 4 J 0 "J 0 (1 + x)(l + y)(l + z)(l + w)(l - xyzw) ' 

(A.2) 


This integral can be done by first making the consecutive and systematic change of 
variables: {x 1 = xyzw, y = y, z = z,w = w}, {x' = x, y' = yzw, z = z,w = w} and 
{x = x,y' = y,z = zw, w = w} 


(n — 1) r°° I" 00 log (x)dxdydzdw 

16vr 4 J 0 J 0 (y + x)(z + y)(w + z)(l + w)(l - x) ' 

(A.3) 


Doing the w integral first we find 

(n — 1) r°° f°° dx logx f°° dzlogz 

16vr 4 J 0 V J 0 (1 -x)(x + y) J 0 (1 - z)(z + y) 

1) j f vr 2 + (log?/) 2 \ 2 _ (n- 1) 

lbvr 4 J 0 y \ 2(1 + y) ) 30 • 


(A.4) 
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B. Calculation of the double sum term 


We are interested in evaluating the term from (j4.12|), which is 


n— 1 


3~ 1 


^(i)*,(i)) 2 ;r«y(i)'Mi)> 2 

3 = 2 j '=1 

1 " 1 / ‘°° dxdydzdw{l — (zw) j_1 )(xt/) j_1 (1 — x)(l — y){l — z)(l — w) 


E 


167r 4 ^ J o (1 — zw)( 1 + x)(l + t/)(l + z){l + iu)(l — x n ){l — y n ){ 1 — z n )( 1 — w n ) 

(B.l) 


One can add to this expression the term j = 1 of the sum, since that will not 
contribute. Therefore carrying out the second sum and taking the n —» 1 limit we 
get: 


(n — 1) f 00 I" 00 dxdydzdw log(xy) 

16tt 4 J 0 "Jo {1 - xy)(l - zw)(1 + x)(1 + y)(l + z)(1 + w) 




dxdydzdw log (xyzw) 

(1 — xyzw)( 1 — zw)(l + a;)(l + y)( 1 + z)(l + w) 


(B.2) 


We observe that these integrals diverge individually for zw —> 1, but not when 
combined. In order to cancel the spurious pole in zw = 1 it is convenient to use the 
same change of variables used in the single sum term 7 : 


fo- 1) 

167T 4 


dxdydzdw log(x) 


/o (1 - tK 1 - z)iv + x)(l + y){w + z)(l + w) 
v dxdydzdw log(a;) 


Jo "'Jo (l ~ x)(l - z)(y + x)(z + y)(w + z)(l + w) 

When added together, the term: cancels the spurious pole. 


(B.3) 


(n — 1) r°° r°° dxdydzdw log(x) 

16tt 4 Jo "Jo (1 - x){y + x){ 1 + y){z + y)(w + z)(l + w) 

Again doing the w integral first, we get: 


(B.4) 


[n 


1 ) 


167T 4 


dz log z 


dx log x 


dy r _ _ 

l + z/io {^-z){y + z)] Q (1 -x){y + x) 


(n 


1 ) 


167T 4 


dy 

1 + y 


7T 2 + {logyf 
2(1 + 7/) 


n 


1) 


60 


(B.5) 


7 In the first integral {x 1 = xy,y = yz' = zw,w = w}, while in the second (a/ = xyzw,y' = 
yzw , z' = zw , w = w} 































C. Calculation of the triple sum term 


The triple sum term given by ([4.18|) 


is 


/ d /? \^ n ' i 4 k 1 

4n ( ^ E E E [<^(1)^(D) 2 (0!(1)0o( 1)) 2 + <0fc(l)0i(l)) 2 (0j(l)0o(l)) 2 

' ' j=3 fc=2 i=l 

+(^(l)0/(l)) 2 (^(l)0o(l)) 2 

+2(0o(l)0j(l))(0i(l)0K 1 ))(0fc( 1 )^( 1 ))(0fc( 1 )0o(l)) 
+2(0o(l)0fc(l))(0fc(l)0i(l))(0i(l)0z(l))(0K 1 )0o(l)) 
+2(0 O (l)0z(l))(0fc(l)0z(l))(0i(l)0fc(l))(0i(l)0o(l)) • 

(C.l) 


Writing all propagators in terms of integrals leads to 


4n 


/ RaRb \ 


(1 — x)(l — y)( 1 — z)(l — w)dxdydzdw 


I o (1 + x)(l + 1 /) (1 + z)(l + w)( 1 — x n )(l — y n )(l — z n )(l — w n ) 


\ r 2 J 167r 4 Jo 

n—1 j —1 fc—1 

EEE [( a; 2/) / ” 1 (- w ) J " fc " 1 + (3w) fc_ ' _1 (xy) J_1 + 

j= 3 fc=2 Z=1 


fc-1 


+2x j - 1 y j - l - 1 z k - l - 1 w k ~ l + 2x K - l y>- K - 1 z 3 - l - 1 w 1 - 1 + 2x l ~ 1 y li ' ^z J " ^ 


r k - 1 od- k - 1 ^j- l ~ 1 nnl - 1 


J-1 ,J-1 


(C.2) 


The procedure to carry out all the sums consists in using the geometric formulas for 
the partial sums, complete the terms to get a sum which starts at k — 1 and j — 1 
respectively, and again apply the geometric formulas for the sums. The result is 

4 „f»V 


(1 — x)(l — 2 /) (1 — 77 ) (1 — w)dxdydzdw 


\ r 2 J 167T 4 Jo Jo (1 + x)(l + y)(l + z)(l + w)(l — x n )(l — y n )(l — z n )(l — w n ) 


71—1 


r 1 

1 — (zw)- 7 1 

(zw)- 7 1 — (xy)- 7 1 

11 — xy 

1 — zw 

zw — xy 


1 — zw 
1 

1 — zw 
2 

l-yz 
2z 


zw 




(. xyzw) 


7-1 


[zw 


it- 1 


1 — xy 

(. xy ) J 1 — ( xyzwY 1 


zw- 


(xr/) 


7-1 


zw — xy 
U - i)0d/) J_1 


V 


1 — xy 

(xy) J_1 — (xw)- 7-1 (xy)- 7-1 — (xyzw)- 7-1 


+ 


z — w 

2 


y — w 1 — zw 

(; yzy~ i — (xwy- 1 (y^y- 1 — (xzy- 1 


+ 


x - 7/ 


wz 


yz — xw 
F -1 — (wx 


|7 —1 


WZ 


yz — xz 
d-i 


(wy) 


7-1 


z — X 


z - y 


(C.3) 
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After careful inspection, it becomes evident that by combining the previous terms 
appropriately, the full integral does not have poles. We can proceed systematically to 
make this fact explicit in order to carry out all the integrals. We found it convenient 
in that sense to combine the first three lines of (|C.3|) and rewrite it as 


(1 — x)(l — y)( 1 — z)(l — w)dxdydzdw 


\ r 2 J 167T 4 


n -1 

E 

3 = 1 


(ZW 


\3~ 1 


ZW 


o (1 + x)(l + y)( 1 + z)(l + w)(l — x n )(l — y n )( 1 — z n )( 1 — w n ) 

(■ %y) j 1 ' 


(xy) 


j-1 


1 — zw 
2 

l-yz 
2z 


3~ 1 


+ 


1 — (zw) 

1 — zw 
(, xy ) j-1 — (xw)^ 1 
V y-w 

(; yz)i~ l — (iw)^ 1 


[zw 


p- 1 - 


xy) 


xy 

(J - 1) . 

zw — xy 

(; xyy^ 1 — ( xyzw)i~ l 
1 — zw 

(: yz ) J_1 — (xz) j ~ x 


j~ 1 


z — w 

2 

x-y [ 


TO 


2/2 

p-i 




(wx) 


i-i 


wz 


2/2 - 
tf- 1 


xz 
- ( wy) j ~ l 


x 


y 


Now take the n —> 1 limit of the full expression: 

2 


-4n(n-l)(^) ^ 


dxdydzdw 


'o 


/ 0 (l + x)(l + ?/)(l + z)(l + w;) 


(1 — zw)(l — xy) 


1 + 


log(xy) log(z;w) log (xyzw) \ 


+ , 

xy 1 


(1 — zw) 

1 


zw 


xy 


log (xy) \og(xyzw)\ 
1 — xy 1 — xyzw J 
\og(zw) log (xy) 


+ 


- zw 
xy 


1 — xyzw J 
d (log (xy) 


1 — zw d(xy) 


xy 


zw 


xy 


2 

y / log(xy) 

log(xta) \ 

1 

flog (xy) 

log (xyzw)\ 

l-yz 

y — w ^ 1 — xy 

1 — xw J 

1 — zw 

\l-xy 

1 — xyzw J 


2z 

\ 1 

f\og(yz) 

\og(xw) ^ 

1 

( log(2/^) 

log(x^) ^ 

w — z 

yz — wx 

11 -yz 

1 — xw ^ 

yz — xz 

11 -yz 

1 — xz ) 


x-y 


log(iy^) log(xu>) 


log(tc^) log (wy) 

z — x \ 1 — wz 1 — xw J z — y \ 1 — wz 1 — wy 


(C.4) 


(C.5) 


We now explain how to compute all the above integrals. 

C.l Integral evaluation 

To carry out the multi-variable integrals of ( |C.5|) it is convenient to analyze the 
integrals in appropriate pole-free combinations: Let’s consider the first three lines of 
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integrals together: 


1 roc roc dxdydzdw 

16vr 4 J 0 " Jo (1 + x)(l + y)(l + z)(l +w) 

log (xyzw)\ 

1 — xyzw J 

1 / log(xy) log {xyzw) \ xy d f log (xy) 

(1 — zw) 2 \ 1 — xy 1 — xyzw J 1 — zw d(xy) \ 1 — xy 

1 / log(^) _ log(xy) \ 1 

zw — xy \ 1 — zw 1 — xy J J 


(1 — zw)( 1 — xy) 


x log (xy) log(zw) 


xy 


1 


zw 


(C.6) 


It is clear from the structure of the integrand that we can performed the next change 
of variables: (x' = xy, y' — y,z' — zy, w' = w) followed by the relabel (x — x',y — 
y',z = z',w = w') of all the integrals simultaneously. That leads to. 


j_ r r ixiz r _*_ r _^_ 

16vr 4 Jo J 0 J 0 (y + x)(l + y)J 0 (w + z)(l + w) 
f_ 1 L log^ log^ _ log(^) \ 

l (1 — z)(l — x) \ l — x 1 — z 1 — xz ) 

1 / logz log(xz)\ x d ( logx\ 

(1 — £) 2 \ 1 — z 1 — xz ) + 1 — z dx\l — x) 

1 { logz _ log(g) \ 1 

z — x \1 — z 1 — X ) J 


We evaluate the y and w integrals, and end up with a finite and well defined double 
integral in x and z 


167T 4 


dx 


logx 
l — x 


dz 


log z 

1-z 


z — X 


(1 — ^)(1 — x ) 


I logx ^ log z log(x^) 


X 


(1 


log z log(xz) 


xz 


+ 


z 

X 


d 


1 — z dx V 1 


xz 
logx 


X 


log£ log(x)\ 
1 — z 1 — X J 


1 

30' 


(C.8) 


This integral can be reduced to a single integral in terms of polynomials of {x, log 
functions as well as poly-logarithmic functions }. The exact value was derived from a 
numerical evaluation. The remaining three lines of integrals in (CJ.5) can be evaluated 
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in a more systematic way by removing the spurious poles one by one, in the integrals 


dxdydzdw 


167T 4 


/ o {l+x)(l + y)(l+z)(l + w) 


2 

J_ 

y ( log (xy) 

log(xw) \ 

i 1 

(log (xy) 

log (xyzw)\ 

- 

i— 1 

1 

o 

y — w V 1 ~ xy 

1 — xw ) 

1 — zw 

V 1 - xy 

1 — xyzw J 


+ 


x-y 


\og(wz) log(xta) 


z — x \ 1 — wz 


1 — xw 


log(wz) log(iuy) 
z — y \ 1 — wz 1 — wy 


2z 

1 

\ 1 1 

flog {yz) 

log(xu') \ 

| 1 , 

flog {yz) 

\og(xz) \ 

r 

w — Z 

yz — wx 

1 1 -yz 

1 — xw ) 

yz — xz 

1 1 -yz 

1 — xz J 


(C.9) 


We consider this line by line, and do a different change of variables in each individual 
integral inside a given line. For example, in 
1 /'°° /'°° dxdydzdw 

16vr 4 J 0 " Jo (1 + x)(l + y)(l + z)(l + w) 

2 y f l°g ( x v) log(xiu) \ 1 / log(xy) log (xyzw) 

1 — yz y — w \ 1 — xy 1 — xw J 1 — zw \ 1 — xy 1 — xyzw 

(C.10) 

We set (V = xy,y ' = y, ...) in the first and third integrals, (V = xw,w' = w, ...) in 
the second, and (x ! = xyzw,y’ = y, ...) in the fourth one. After removing the primes 
we have 



1 f°° f°° dydzdw f°° (log x \ 

Te^Jo '"Jo (l + y)(l + z)(l + w)J 0 dX {—x) 


2 

\ V 1 

( 1 

1 ^ 

1 

( 1 

1 

l -yz 

_y - w 

\y + x 

W + X J 

1 — zw 

\y + x 

yzw + x J 


(C.ll) 


with these simple steps we cancel all the poles explicitly by simply adding the different 
fractions 

1 I* 00 f°° dydzdw f°° /loga:\ yw 

8tt 4 Jo Jo (1 + I/)! 1 + -)(! + w) Jo ^ V 1 - x ) iv + x )(w + x){yzw + x) ' 

(C.12) 


As a further simplification we do the next change of variables: (V = zyw,y' = 
yw,x' — x, w' — w), and carry out the integrals on w and z respectively 


1 r° / log x \ dwyw r 

8tt 4 Jo X \1-x)J 0 V J 0 (1 + w)(w + y)(y + xw)(w + x) J 0 

1 f°° f°° dxdyy /loga;\ /logo: — log/ log(x 2 ) 

8vr 4 Jo Jo (x-y)(x 2 - y) \1 - x) \ x-y J \ 1 - x 

1(1 15 \ _ ( 1 1 \ 

~ _ 45 V2 + 4^2 J “ _ V90 + 12^2 ) ■ 


dz 


(y + z)(z + x) 
log y 


(x + y) 


1 -y 


(C.13) 
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After following the same steps as the previous evaluation, we get for the final integral 
expressions for the second and third lines of QC.9|) : 


1 f°° f°° dxdy / log x \ / logo: 

8vr 4 J 0 J 0 (x-y) 2 \l-x) \l-x 
in 15 \ _ /1 i \ 

45V2 + 47T 2 y V 90 + 127r2 / 


log y\ ( logx log y\ 
l-y)\ X l-x y i-y) 

(C.14) 


and 


87T 4 


dw 


login 
1 — w 


“~45 ,1 + 


15 

2^2 


respectively. The hnal value of the integral in 


is: 




(C.15) 


“(^ + A + 3?)- (CU6) 

With this result at hand we can write the final result for the triple sum term, which 
is 


71—1 

S C fj'lO C fj'lO = 4n ( n 



12 1 \ 
30 + 45 + 3^V ’ 


(C.17) 
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